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Abstract. We use energy minimization principles to predict the structure of a 
decagonal quasicrystal - d(AlCoNi) - in the Cobalt-rich phase. Monte Carlo 
methods are then used to explore configurations while relaxation and molecular 
dynamics are used to obtain a more realistic structure once a low energy 
configuration has been found. We find five-fold symmetric decagons 12. 8A in 
diameter as the characteristic formation of this composition, along with smaller 
pseudo-five-fold symmetric clusters filling the spaces between the decagons. We 
use our method to make comparisons with a recent experimental approximant 
structure model from Sugiyama et al (2002). 



PACS numbers: 61.44.Br, 61.50.Lt, 61.66.Dk, 64.60.Cn 

1. Introduction 

This paper reports structural predictions for the decagonal quasicrystal d(AlCoNi) in 
the Cobalt-rich ('basic Co') phase , of approximate composition Al7oCo2oNiio and 
atomic density 0.068 atoms/A 3 , using simulations to minimize energy, using the same 
methods (and codes) as previous papers by Mihalkovic, Widom, and Henley [1, 2, 3] 
on the 'Ni-rich' phase ; this constitutes a first test of the transferability of that 
approach to other compositions that are described by other tiling geometries. 

The Aluminum-Nickel-Cobalt alloy d(AlCoNi) is the best studied and perhaps 
the highest quality of equilibrium decagonal quasicrystals. Decagonal d(AlCoNi) 
has about eight modifications (which we call 'phases ' here), each existing in a 
tiny domain of the phase diagram. [4, 7, 5, 6, 8]. The compositions 'basic Ni' 
(around Al7oCoioNi2o), and 'basic Co' are the best-studied phases , having the 
simplest diffraction patterns. (The other phases show superstructure peaks, indicating 
modulations.) Several crystal approximants related to 'basic Co' are known [9], one 
of which has a solved structure [10]. 

We find that (in contrast to the 'basic Ni' phase ), the framework of 'basic Co' 
at T = is a network of edge-sharing 12. 8A diameter decagons (placed like the 
large atoms in a "binary tiling" quasicrystal). A strong (but not always simple) Co/Ni 
ordering is found. Our model reproduces most, but not all, of the atomic positions 
in W(AlCoNi). (A complementary brief account of our results is in Ref. [11], and a 
more complete account is in preparation [12].) 
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2. Methods 

Our main assumption is that the d(AlCoNi) structure is well approximated by a 
stacking of equally spaced two-dimensional tilings built from the Penrose rhombi 
under periodic boundary conditions. The chief experimentally determined inputs 
were the values c =4.08A for the stacking period and a =2.45A for quasilattice 
constant (edge length of small rhombi). It should be noted that the actual period 
for the Co-rich phase is generally believed to be 2c; we used c because (i) it is the 
simplest starting point, and was used in the prior study of the Ni-rich phase [1,2] 
(ii) so long as the atoms are fixed on discrete sites (see Section 3), it turns out that 
period c is optimal even when a cell with more layers is allowed; (iii) even in our 
final, relaxed period 2c structure (Section 4), most atoms do repeat with period c. 

This framework of decorated tilings takes advantage of the fact that all the 
known atomic structures of decagonals are Penrose-tiling-like. Further, even in 
random tilings Penrose rhombi admit "inflation" constructions that can relate a tiling 
to another one using rhombi with edges enlarged by r = (1 + /2 w 1.618 (the 
golden ratio), which gives a convenient way to build a chain of connections from the 
atomic level to large tiles so as to describe d(AlNiCo) on large scales. Fig. 1 includes 
two ways of subdividing a unit cell into a tiling of rhombi (with edges ao and rao, 
respectively). 

We also assume atomic pair potentials between Al, Co, and Ni as derived ab- 
initio using Generalised Pseudopotential Theory [13] (GPT), but modified using 
results from ab initio calculations to approximate the sums of the omitted many-body 
interactions when transition metals (TM) are nearest neighbours. [Fig. 1 of Ref. [1] 
plots these potentials.] The potentials are cut off by a smooth truncation past 7.0A. 

We know [1, 2] that the strongest interaction is the Al-TM first potential well. 
The number of Al-TM neighbours is maximised and, as a corollary, that of TM- 
TM nearest neighbours is minimised (but TM-TM second neighbours, at ~ 4A 
separation, are common). However, this scarcely constrains the structure since the 
weak Al-Al potential allows enormous freedom in placing Al. In practice, one finds 
in Al-TM quasicrystals that a rather rigid TM-TM network forms, with separations 
near to the strong second minimum of the TM-TM pair potential; the Al atoms fill 
the TM interstices in a more variable fashion. Subtle details of the density and 
composition, as well as the small differences between Co and Ni in the pair potentials, 
decide which structure optimises the energy. We understand only a few of these 
factors in detail for the 'basic Co' phase, so in this paper we will only describe the 
atomic structures we found and do not attempt a microscopic rationalization of them. 

Our program uses Metropolis Monte Carlo (MMC) simulation to perform atom- 
atom and atom-vacancy swaps within a site list placed on the Penrose tiles. The first 
stage simulations include an additional degree of freedom in the form of 'tile flips.' 
These 'tile flips' obey the MMC algorithm probabilities and conserve tile number, 
atomic species number, and the collective outline of tiles rearranged. 

Our procedure is to anneal from a high temperature j3 ~ 4-10 eV -1 to a low 
temperature at which most degrees of freedom are frozen out, (3 « 20 eV _1 , using 
a step size 6(3 ~ 0.5-1 eV _1 . (Here j3 = (fc B T) _1 .) At each temperature step, 
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Figure 1. The 32x23x4 unit cell. The left half is subdivided into ao scale rhombi 
(in two ways in part of the cell [dotted lines]). The right half of the cell is subdivided 
into tcio scale rhombi. The atomic configuration shown was obtained from the ran 
scale simulation. (The atomic locations are also valid for the ao scale decoration, 
but configurations from ao scale simulations almost always have more defective 
occupations. )Insets at top left and top right display the candidate sites on the ao 
and rao scale tiles, respectively. Filled large circles depict atoms in one layer and 
filled small circles are those in the other layer. Black denotes Nickel, dark gray 
Cobalt, and light gray Aluminium. Unoccupied sites of the rao scale tiling are 
shown by very small empty circles; those of the ao scale tiling are not shown. 
Four 13A Decagon motifs (described in Sec. 3.1) are present: e.g. one is seen 
decomposed into ten rao edge rhombi near the cell's upper right corner. 

approximately 2000 atomic swaps per site are attempted. 
2.1. Multiscale Procedure 

The strategy of our series of simulations loosely copies that reported previously [1,2] 
on the Ni-rich phase of rf(AlNiCo), and can be thought of as a sort of multiscale 
modelling. It proceeds in the following stages: 

1. Begin with Monte Carlo simulations allowing atom swaps and 'tile-flips' 
with a tiling of relatively small Penrose rhombi; note prominent features and 
formations that recur in all the low energy structures from a sufficient number 
of independent runs. 

2. Take the observations from the small-scale tiling and promote them to rules: 
i.e. make the formations inferred in the first stage into fundamental objects of a 
larger scale tiling. Unused (or underused) degrees of freedom are removed. 

3. Perform new Monte Carlo simulations on the larger-scale tiling to search for 
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new low energy structures. This larger- scale tiling is generally more efficient at 
finding lower energies because of the reduced degrees of freedom. 

4. Verify that the degrees of freedom removed in stage (2) were indeed unnecessary 
by restoring some of them, judiciously. 

5. Repeat for larger and larger scale tilings. An eventual goal is to extract a 'tile 
Hamiltonian' [14, 15] which is the bridge to modelling long-range order and 
diffuse scattering. That means ascribing the atom-atom energies to an effective 
interaction among neighbouring tiles, so that tiles are the only remaining degree 
of freedom. 

It is also necessary to perform relaxations/molecular dynamics to escape the 
biases which may be introduced by the discrete site list used up through stage (4). 

3. Results 

3.1. Small-Scale tiling (edge a — 2A5A) 

The first stage of our tiling simulations uses small Penrose rhombi that have ao = 
2.45A edges. The unit cell simulations had dimensions 32.0lA x23.25Ax4.08A, 
which we shall call the 32x23x4 cell, shown in Fig. 1. The structure is taken to 
repeat after two layers of the Penrose tiling. At 0.068 atoms/A 3 , there are about 200 
atoms in the unit cell. 

Our low energy configuration generated at this level showed distinct rings of 
5 TM and 5 Al atoms surrounding an Al atom and surrounded by an 10 Al atom 
decagon with an edge length of ao. These 21 -atom motifs were evidently energy 
favourable. Linear regressions of the count of these 21 -atom clusters versus energy 
revealed a noisy, but consistent correlation of w — 1.0 eV per 21 -atom cluster. 

However, the ideal structure did not simply maximsze the number of 21- 
atom clusters: though they may be placed adjacent to pack six clusters per cell, 
that was typically not observed in the low-energy configurations. Instead, further 
MMC annealing (starting from a low-energy configuration) found that a larger 
cluster emerged as the characteristic motif at this composition: a decagon with an 
edge length of too w 4.0A (hence a diameter of ~ 12. 8A), which we will call 
' 13A Decagon' (13AD). This object contains the 21-atom motif at its centre, ten TM 
atoms at its vertices, and Al atoms irregularly interspersed along and within its edges. 

We use the following nomenclature for the components of 13AD, from its centre 
outwards. The centre of the 13AD is an Al atom. The 5 TM 5 Al ring is known as 
the 'first ring.' The TM atoms in the first ring are in the same layer as the central 
Al atom. The 10 Al ring is known as the 'second ring.' The 10 TM decagon and 
the Al atoms along and within the edges are collectively known as the 'third ring.' 
Miscellaneous Al sites occur between the second and third rings, which we designate 
collectively as the '2.5 ring.' 

In the 13 AD, the first ring TM sites are mainly occupied by Co and the third 
ring (decagon vertex) TM sites also tend to Co occupation, but there is a very strong 
dependence on the local environment that distinctly favour Ni in certain sites. [This 
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became clear at the second stage of simulation, Subsec. 3.2.] The Al positions of 
third-ring and 2.5 ring Al atoms are closely correlated and the rules were not resolved 
from the 2.45A tiling simulation. 

A second type of motif fills the interstices between the 13ADs: it consists of a 
10-atom ring and an Al at the centre, much like ring 1 in the 13AD cluster, except 
this kind of motif has mixed Al-TM occupations for the 5 candidate TM sites. In 
general, about 3 sites out of the 5 candidates are filled with TM atoms, preferably 
Ni. We name these 1 1-atom clusters "Star clusters"; a few of these are contained in 
Fig. 1. 

Perfect examples of either cluster were rarely observed in the first (2.45 A) level 
of simulation: their decorations were induced by a 'consensus' or average over many 
defective examples. Thus, the rules for Al/TM occupation on the Star cluster TM 
sites, and for Ni/Co occupation on all TM sites, were still incompletely known from 
this stage. 

3.2. Inflated Tiling (V 3 ao — * Tdo) 

In keeping with our 'multiscale' methodology, we reassessed the degrees of freedom 
necessary for the next level of simulation. The dominant geometric object is a 
decagon with edge rao ~ 4A, but we found it economical to represent its decoration 
by candidate sites using a thin 4A rhombus and the second fat rhombus shown in the 
right inset of Fig. 1; similarly, the Star cluster decoration is represented using the 
first fat rhombus. 

We assumed that the minimum-energy structure has a maximum density of 
13 AD clusters. A Monte Carlo simulation was carried out (admitting only tile flips) 
on the ensemble of rao rhombus tilings using an artificial tile Hamiltonian favouring 
the creation of nonoverlapping 'star decagons' (of edge rao) containing a fivefold 
star of fat rhombi. % It was found that the ground state configuration was always 
a 'binary tiling' [16, 17] of rhombi with edge length r 3 ao =10.4A with the 'large 
atom' and 'small atom' vertices replaced by 13ADs and Star clusters, respectively. 
(The constraints on possible decagon separations imposed by building them from rao 
tiles are consistent, in fact, with those due to the potentials [12].) 

We designate a fixed site list on the rao rhombi, as shown in the right half 
of Fig. 1, designed to match the most frequently occupied sites observed in the 
low-energy states from out a scale simulations. This tiling, unlike the a tiling, 
consists of only one layer of Penrose tiles, each of which has site decorations in two 
layers spaced 2.04A in the c direction. Note how the two versions of the fat Penrose 
rhombus have different site lists. We use the 32x23x4 cell, as before so we may 
compare the energies for systems with exactly the same density and composition. 

At this second stage, our MMC runs had purely lattice gas moves on a fixed 
tiling (no tile flips). These were able to find lower energy configurations at a much 
faster rate than the first-stage (a scale) simulations at the same composition and 
atomic density, due to the much decreased site list. Hence, the occupations of many 
types of sites were discerned with greater accuracy, such as the TM sites in the 

| Star decagaons were similarly maximized in Ref. [18], except they could overlap in that case. 
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13A Decagon (already described in Subsec. 3.1). Furthermore the tendencies of 
more variable sites (Al in rings 2.5 and 3, and TM in the Star cluster), which will 
be described in a longer paper [12], started to come into focus. Finally, despite 
having fewer allowed configurations, the rao tiling simulation consistently found 
lower energies, which serves as our post hoc justification for assuming 13AD clusters 
in the binary tiling geometry and the reduced site list. § 

In the prior case of the 'basic Ni' phase , the decoration model [1] was 
implemented as a deterministic tiling. A marked contrast in the present 'basic Co' 
case is that we cannot impose a simple rule that fixes the chemical occupation of 
each site. On certain TM sites, it is difficult to resolve the occupation (Co/Ni); 
certain Al sites are also variably occupied (Al/vacant) depending on the environment. 
Presumably, with a sufficiently thorough understanding of our model, we could 
formulate a deterministic decoration rule on the binary tiling. It is possible, however, 
that the energy differences among some competing structures are too small to be 
visible in a reasonable simulation, or to influence the real properties at any accessible 
temperature. 

3.3. Cluster orientation 

The physical accessible candidate sites in the rao decagons (with our decoration) 
in fact have a 10 point group symmetry, which is not broken by the binary tiling 
geometry. The fivefold symmetric 13AD cluster breaks this symmetry, by the layer 
on which the central Al sits and the Al/TM alternation on the first ring. Thus, a 
major obstacle to writing a deterministic decoration is that the relative orientations of 
13ADs must be specified, which depends on subtle interaction energies between 
them. Experimentally, <i(AlCoNi) with our Co-rich composition was observed 
to order in a fivefold symmetry, wherein all the Decagons are oriented the same 
way [19, 20]. 

As a test, we compared the lowest-energy configurations resulting from MMC 
simulations on identical ra tilings in which neighbouring 13ADs were forced 
either to have always identical or always opposite orientations. [This is controlled 
by the orientation of the decagon of ra rhombi in the tiling that gets decorated.] we 
evaluated the minimum-energy configurations from ~ 60 independent runs of each 
type. We used the 32x23x4 unit cell, using a fixed composition Al7 Co 2 oNii , but 
repeating the tests for a a series of atom number densities n. 

The energy difference between these two orientations was density dependent. 
At n w 0.068 atoms/ A 3 the two orientational patterns are practically degenerate 
(though with visibly distinct atom configurations.) From n =0.068 A~ 3 to 0.074 A -3 
- a range which includes the most realistic compositions - the pattern with uniform 
orientations has the lower energy, by up to 4 x 10~ 3 eV/atom, but this difference 
disappears again at n «0.074 A~ 3 . 

We find this energy difference arises from different ways of filling the possible 
TM sites in the respective orientation schemes. In the uniform-orientation scheme, as 

§ In additional tests, certain candidate sites that are absent in the inflated tiling were restored 
systematically to the rao tiles, and it was verified this had a negligible effect. 
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the TM density is raised (as part of the total density), TM atoms fill the ring 2.5 in the 
13A Decagon before they exceed ~ | filling in the Star cluster TM sites; whereas in 
the alternating-orientation case, the ring 2.5 is filled to a lesser degree with TM while 
the Star cluster TM sites become overfilled. [In Ref. [11], we gave an explanation 
how these effects could favour alternating orientations.] But the results just presented 
are valid only for the (unphysical!) case of atoms confined to hop on fixed ideal sites. 
When displacements are allowed, so that many atoms "pucker" away from the flat 
layers (see Sec. 4), then the uniform arrangement is preferred more robustly. 

4. Beyond the Discrete Site List 

Up to this point, we reported simulations using fixed atom site positions. We 
can make our configurations more realistic by subjecting them to relaxation and 
molecular dynamics (MD). The removal of the tiling and discrete site list leads to 
effects such as out-of-layer "puckering" and so-called 'period doubling,' whereby 
certain atoms relax into positions that violate the c period of a single unit cell, but 
instead are periodic with respect to a unit cell with c' = 2c = 8. 16 A. Exactly such 
distortions are familiar in the structures of decagonal phases and approximants. 

Our standard cycle for such 'off-site' studies is a relaxation to a local energy 
minimum (i.e. to K), followed by an MD cycle beginning at ^600 K and ending at 
^50K in increments of 50K. The MD results are then relaxed again to OK. This 
protocol will be denoted Relaxation-MD-Relaxation (RMR). Upon relaxation to 
T = 0, we find that the TM atoms are quite immobile and move only slightly. The 
Al atoms, however are subject to displacements as large as 1.5A. After molecular 
dynamics and re-relaxation to T = 0, we find a few Al atoms be further displaced, 
but to similar sites so that no systematic difference is apparent in the overall pattern. 

4.1. Comparison with Experimental W-AlCoNi 

Major diffraction-based structure models were available for 'basic Ni' d(AlNiCo), to 
which the simulation predictions could be compared, but no such model exists for the 
'basic Co' d(AlCoNi). However, the structure of the 'W phase' crystal approximant 
of 'basic Co' has been determined by Sugiyama et al [ 1 0] , and we may apply what we 
have learned so far - in particular the period doubling and puckering - to predict its 
structure. [This study has inspired a more detailed modelling of W(AlCoNi), based 
on ab-initio energies rather than pair potentials [21].] 

Using exactly the same tiles and interlayer spacing as in our other simulations, 
our unit cell is 23.25Ax39.5606Ax8.158A, which differs by less than 1% from 
the experimentally determined lattice constants. The approximate atom content 
implied by Sugiyama's structure solution and Co:Ni ratio is Al3 85 Coii3Ni 38 , 
with uncertainties arising from the partial and mixed occupations. Adopting this 
composition as our ideal, we have Al71.sC021.1Ni7.! at a density 0.0714 A~ 3 . 

We applied our mock Hamiltonian from Sec. 3.2 to generate a tiling consistent 
with the r 3 ao binary rhombus tiling. The optimum configuration is unique (modulo 
symmetries) and has four 13ADs, in an arrangement which turns out to be the same 
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as observed in in W(AlCoNi). Cluster orientation comparisons were like performed 
like those of Subsec. 3.3, but using the relaxed energies (without MD). We found 
that, between the configurations with alternating and uniform cluster orientations, 
the latter had a lower relaxed energy, in accordance with the experimental W-AlCoNi 
structure. The atomic locations generated on this unit cell using the 4A rhombus 
decoration site list is approximately consistent with the experimentally determined 
sites. 

Our MMC simulations are able to capture the gross features of W-AlCoNi with 
excellent accuracy. Since this approximant is not far from the 'basic Co' decagonal 
composition, we are not surprised to find 13ADs and Star clusters as the major 
motifs, arranged in a manner consistent with our binary and rao tilings, consistent 
with the observed W-phase. The differences in the exact atomic locations lie mainly 
within the highly context dependent 2.5th and 3rd rings. We find that the increase 
in density from our original simulations results in a more dense occupation of these 
rings. Atomic configurations found after RMR are in even better agreement. Fig. 2 
compares one layer from experimental data to our RMR results; the match in the 
other layers are equally good. 

However, a significant number of Al atoms in our RMR picture are in 
disagreement with the experimental refinement, because certain puckered atoms 
are not present in our simulations. The top right of Fig. 2 contains a pentagon of 
displaced Al atoms centred upon the missing Al atom. These extremely puckered 
Al atoms are perhaps located in energy minima too far away for our MD program to 
traverse. The existence of such atoms is one of the greatest flaws for beginning with 
a discrete site simulation. 

5. Conclusion 

We find that in the approximate composition (AlyoNiioCo2o, the structure of 
d(AlNiCo) is dominated by the formation of the 13AD clusters and the Star 
clusters which complement them. We find that the gross features are robust under 
variations of the composition and densities of ~ 4%; these will strongly affect a few 
details of the sites and we have not established the correct answers for those atoms. 
Our method of tile decorations, and successive reductions of the degrees of freedom, 
affords an enormous speedup compared to brute-force molecular dynamics, which 
would get stuck in glassy configurations, at the price of a few plausible assumptions. 

By applying our results to the approximant W-AlCoNi, we demonstrated both 
the validity and the pitfalls of our approach. The pitfall lies in that fact that our 
explorations are conducted using a fixed site list in two atomic layers, whereas the 
real structure has four layers and some atoms are strongly puckered out of them. 
Even though most puckered atoms may be represented as relaxations from ideal 
sites, the optimal puckered structure might be derived from a fixed-site structure of 
comparatively high energy that our simulations would pass over. Reliable answers 
can not be obtained by pure numerical exploration, but seem to demand some 
physical understanding of the puckering (and other displacements), as we have begun 
to do in Ref. [11]. The context dependences evident in the 'basic Co' show that 
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Figure 2. a) The top figure is one of the mirror layers of W-AlCoNi. b) The bottom 
figure is the puckered layer of W-AlCoNi. The empty circles are experimentally 
determined atom sites, the filled circles are those extracted from our simulations. 
Dark gray filled circles are TM, light gray filled circles are Al. Empty circles 
with 'teeth' are TM, smooth empty circles are Al. The large empty circles are 
experimental atomic sites with mixed Al-TM occupancies. The size of the simulated 
(filled) atoms may vary due to the visual scheme depicting distance in the c-axis 
dimension. 



beyond a point, the errors in defining the pair potentials will surely exceed the energy 
differences due to some swaps of species, which is a second pitfall that may be 
alleviated only by checking against experimental phase diagrams. 
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